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• Abstract: The Yang-Mills gradient flow is considered on the four dimensional torus for 

SU{N) gauge theory coupled to Nj flavors of massless fermions in arbitrary representations. 
The small volume dynamics is dominated by the constant gauge fields. The expectation 
value of the field strength tensor squared TrF^j^F^j^(t) is calculated for positive flow time 
t by treating the non-zero gauge modes perturbatively and the zero modes exactly. The 
finite volume correction to the infinite volume result is found to contain both algebraic 
and exponential terms. The leading order result is then used to define a one parameter 
family of running coupling schemes in which the coupling runs with the linear size of the 
box. The new scheme is tested numerically in SU{3) gauge theory coupled to Nj = 4 
flavors of massless fundamental fermions. The calculations are performed at several lattice 
spacings with a controlled continuum extrapolation. The continuum result agrees with the 
perturbative prediction for small renormalized coupling as expected. 
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1 Introduction and summary 

The Yang-Mills gradient flow - or Wilson flow ~ has proved to be a useful tool in lattice 
gauge theory. In the context of the Nicolai map it was studied in [1]; see also [2] for an earlier 
appearance. A systematic investigation, including suggestions for possible applications, has 
appeared relatively recently [3-5]. See also [6]. The first concrete very useful application 
of the flow for high precision setting of the physical scale in QCD simulations has been 
presented in [7]. The flow in QCD applications has so far been considered in infinite volume 
which is most appropriate for low energies. 

In the present work the flow is calculated on the four dimensional torus, i.e. in a 
finite four dimensional box. The motivation for doing so is to obtain a new running 
coupling scheme in which the renormalized coupling runs with the linear size of the box. 
In principle the original infinite volume flow can also be used for defining a renormalized 
running coupling gR{q) with q = where t is the flow time, but the control of finite 

volume corrections is an additional concern in this case. This issue is eliminated if the 
running gR{L) is with the linear size L. In particular, a step scaling analysis can be 
performed [8, 9]. 

Due to asymptotic freedom perturbation theory is reliable for small volumes hence the 
appropriate framework is the small volume expansion or femto world [10-14]; see also [IS- 
IS]. The usual complication associated with calculations in the femto world is the presence 
of gauge zero modes which dominate the dynamics and are not Gaussian. They need to 
be treated exactly while the gauge non-zero modes can be integrated out in perturbation 
theory. As will be shown, the contribution of the non-zero modes renormalizes the bare 
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coupling according to the 1-loop /3-function and generates an effective action for the zero 
modes. 

The quantity which turns out to be the most useful for our purposes is the one that 
has already been calculated in infinite volume in [3], namely the field strength squared at 
t > flow time, 



E(t) = - TTF^,F^,{t) (1.1) 



(see appendix A for our conventions). The expansion of its expectation value in flnite 
volume is our main result and to leading order in the MS scheme it is given by 

ifE{t))=gl{^,f-^^^^{l + 5) (1.2) 

where is the dimensional regularization scale, gjiifJ-) is the renormalized coupling in the 
MS scheme. The correction factor 6 = 6a + 5e is a, sum of algebraic and exponential terms, 

64t27r2 



3L4 

4 = # (exp (-^)) - 1 = 8exp (-|!) + 24exp (-^) + . . . , (1.3) 

and where ^{q) is the standard Jacobi elliptic function (normally called '^^{q), see appendix 
A for details). Indeed, the inflnite volume result in [3] is reproduced. 

Equation (1.2) can be used to deflne a running coupling gR{L) which will run with the 
linear size once the dimensionless combination c = ^/si/L is held fixed and fi = 1/L is set. 
Different choices for c correspond to different schemes. 

The organization of the paper is as follows. In section 2 the small volume expansion 
is given on and the finite effective action for the gauge zero modes is calculated by 
integrating out the non-zero modes to 1-loop. In order for the presentation to be self- 
contained all details are spelled out although the methods are by no means new. In section 
3 the gradient flow is considered and the expectation value of the quantity E{t) is calculated 
by again treating the non-zero modes in 1-loop perturbation theory and using the previously 
obtained effective action for the zero modes. The result is then used in section 4 to deflne a 
renormalization scheme for the gauge coupling. As an illustration of the method, numerical 
simulations are used to compute the running coupling in SU{3) gauge theory coupled to 
Nf = 4 massless quarks in section 5. Finally we close with conclusions and provide an 
outlook in section 6. 



2 Small volume expansion 

On the four dimensional Euclidean torus with periodic boundary conditions for the 
gauge field the zero momentum (constant) gauge mode is separated from the first non-zero 
momentum mode by the gap 2tt/L and dominates the low energy small volume dynamics 
[10]; see also [11-18]. This dynamics is non-linear because of the quartic interaction and 
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needs to be treated exactly while the dynamics of the non-zero modes can be treated 
perturbatively. Correspondingly the gauge field is split 

A^{x) = B^ + Q^{x) , j (fxQ^ix) = (2.1) 

into the zero mode and non-zero modes Q^{x). The action for Nf flavors of massless 
Dirac fermions in representation R is 

S = -^ [ d^xTiF^^F^, + I d^xijflpijf . (2.2) 

where go is the bare coupling constant. The boundary condition for the fermions is assumed 
to be anti-periodic in at least one direction. It is convenient to introduce 5^-|--B^ = D^{B) 
acting in either the adjoint or representation R depending on whether it is applied to a 
gauge field or fermion. 

Gauge fixing is only required for the gauge non-zero modes and a convenient gauge 
choice is the background gauge x = Dfj^[B)Q^ = 0. The constant gauge transformations 
do not need to be fixed as their volume is finite. 

Neglecting interactions which are higher order in and the ghost field one obtains 
the leading order Faddeev-Popov operator as Df/_{B)^ which is understood in the adjoint 
representation and acts on ghosts without zero- modes. The corresponding effective action 
for the zero mode B^ is then 

Sg,,{B) = - In det {D ^{Bf ) . (2.3) 
The quadratic term in Qf^ from the gauge action is 

-^2 / d^xTvQ^ {Dp{Bf5^, - D,,{B)D,{B) + 2[B^, B,]) Q, . (2.4) 

A convenient way of implementing gauge fixing is by adding /'^9o to the action which 
allows integrating out the field without the gauge constraint. The effective action from 
this bosonic integral is then, 

Sq{B) = i Indet {Dp{Bf6^, + 2[B^,B,]) . (2.5) 

In the fermionic action one may neglect the interaction between the fields and the 
fermions. To leading order one obtains the effective action 

Sf{B) = - Indet (0(5))^/ = - Indet iD,,{Bf + -a^,[B^,B,]\ , (2.6) 

where a^^ = [^^,^y\/2. Here the operators act on fermions with the appropriate boundary 
condition. The various determinants will be evaluated using dimensional regularization 
and all subsequent calculations are done in dimension d = 4 — 2e. 
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The total effective action after integrating out the gauge non-zero modes, the ghosts 
and the fermions is then 



Seff{B) = -^^^L_Tt[B„ + Sq{B) + S,h{B) + Sf{B) , (2.7) 

where the first term is the tree level action for the constant mode and is the scale of 
dimensional regularization. 

Now we will proceed to evaluating the various determinants. They will be Taylor- 
expanded in Bfj^ and we will see later that it is enough to expand them to fourth order for 
our purposes. Higher orders in will correspond to higher orders in the renormalized 
coupling. The expansion is around the free B^ = determinants and these (infinite) 
constants are dropped as usual. 

The derivatives in Sq and Sgh are replaced by 2Tiin^/L where are integers and 
7^ 0. In Sp the derivatives are replaced by 27ri(n^ — kfj)/L where A;^ is 1/2 in all anti- 
periodic fermion directions and the rest of its components are zero. We will assume /c^ 7^ 0. 
It is furthermore convenient to introduce the hermitian matrices = LB^/27ri. 

Straightforward calculation yields that up to fourth order in the following holds 

SgiC) + SghiC) = TradlogiD^iCf) +^Trad[Cf,,C,f , (2.8) 
where the traces are in the adjoint representation and 

7 = E;^- (2-9) 

Similarly, the fermionic contribution to the effective action up to fourth order in is 

Sf{C) = -2Nf (ttr logiD^iCf) + 2^TYr[C^, C,f^ , (2.10) 
where all traces are in the representation R and 

n ^ ' 

Equations (2.8) and (2.10) show that only the Laplacian is needed in the background of 
in arbitrary representation and with arbitrary boundary condition in order to evaluate 

the full effective action. 

First, let us evaluate all determinants with periodic boundary condition and get back 

to the case of non-trivial boundary conditions for the fermions later. Explicit calculation 

yields up to fourth order in C^, 

- TVh \og{D,{Cf) = 6 ^ Tr^C^ + 7 ^ Tr^^^' + 

a 2d 

+ 4 ^ !h^!hi^TvRC,C.C,C, , (2.12) 
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where the new constant 6 has been introduced and = C^C^ and = {C^C^ 
50(4) invariant combinations. It is useful to define two more constants a and /3 by 

nj^O n=/=0 n=iO 

Using these the following is easy to show, 
E IV^^JpI^Tr^C'.aCpa = (a - 3/3) J] + /? ('3Tr«C^ + ^Tr^[C^, C.]^) . (2.14) 

Since the torus breaks rotations the S'0(4)-breaking first term on the right hand side is 
allowed. Combining equations (2.12) and (2.14) we obtain, 



- Tr^ \og{D^{Cf) = 6 'jz^T^^rC^ + 4(a - 3/3) J] + (2.15) 

+ (l2/3-7|^) TTRC^ + 2/3Trn[C,„C,]^ . 

In appendix B it is shown that even though a, (3 and 7 are all divergent the combinations 
appearing above for the terms that were not present at tree level, namely C^, and 
^^C^, are all finite. Only the coefficient of [C^,C^]^ is divergent. 

Now the full effective action (2.7) is easily written down using (2.15) in the adjoint 
representation together with (2.8) and in representation R together with (2.10). The 
traces of the product of two Lie algebra elements in different representations can be all 
converted to the fundamental representation using the trace normalization factors T[R) 
via Tr/j(- •) = 2T{R) Tr(- •). Let us first collect the terms proportional to Tr[C^, C^]'^ which 
is the only divergent term. Using T{ad) = N and the poles of /3 and 7 from appendix B 
we obtain, 

Seff{C)\,,^ - \-g^(ilLfe + ^^1*^) Tr[C^, C.] (2.16) 



Clearly, by introducing the renormalized coupling gR{i-i) of the MS scheme, 

1 _ 1 ^N-jT{R)Nf 



(2.17) 



in place of the bare coupling go a finite effective action is obtained. Going from MS to MS 
scheme only modifies the finite terms. 

Up until this point the momentum sums corresponding to the fermions were computed 
with periodic boundary conditions, however we are interested in fermions that are anti- 
periodic in at least one direction. Instead of the coefficients a, (3 and 7 we should have 
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considered a:^{k), /3f^,^{k) and ^{k), 

n ^ ' 

/3.^(fc) = E^^HSw^ (2.18) 

n ^ ' 

where kf^ ^ determines the boundary conditions. However, it is easy to see that the 
differences a^{k) — a, jif_iu{k) — j3 and 7(/c) — 7 are ah finite. This is expected because UV 
divergences are insensitive to boundary conditions. Hence once the UV divergences are 
canceled only the finite terms can be effected by the change of boundary conditions. 

Summarizing this section, a finite effective action is obtained for the gauge zero modes 
of the form, 

SeffiC) = -^^Tr[C7^,a]2 + (2.19) 
+uiTrC72 ^ ^^Tr^C^ + u-iT,adC^ + 

where the finite expressions ui, ... ,7x5 depend on N, Nj, R and the boundary condition 
for the fermions. These are all known although in a bit cumbersome form. Their values 
will not be important for what follows, the only property we need is their finiteness. From 
now on we set fi = 1/L. 



3 Yang-Mills gradient flow on 

Now that a finite action is obtained for the gauge zero modes let us turn to our ob- 
servable of interest, the field strength squared E{t) at positive fiow time (1.1). It will be 
evaluated by treating the gauge non-zero modes in perturbation theory and the zero mode 
exactly, similarly to the effective action. Let us first write down the Yang-Mills gradient 
flow, 

^ = «„F„„. (3.1) 

Using the decomposition (2.1) we obtain a coupled fiow for the zero and non-zero modes. 
After taking into account gauge fixing and dropping terms higher order in we arrive at, 

[B,,[B,,B^]] (3.2) 
{Dp{Bf5^, + 2[B^,B,])Q, . 



dB^ 
dt 

dt 
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Since we are interested in a perturbative expansion let us rescale — gnQ^. The 

1/2 

consistent rescaling of the zero mode is — )■ (7^ B^. After the rescahng the gradient flow 
becomes 



^ = g^[B,^[B,,B,]] (3.3) 
= AQ^ + Oigjl ) . 

Clearly, to leading order in the coupling the zero mode is constant B^{t) = B^ and the 
solution for the non-zero mode is 

Q,,{t) = e*^Q^(0) . (3.4) 

In the path integral one integrates over the fields at t = 0, i.e. Qfj,{0) and i?^. 

The rescaling also effects the observable E{t) and keeping the leading order term only 
we obtain, 

2 2 

m = -^Tv[B^, B,f + ^Tr Q^e^*^ (A^ " 9^,8,) Q, (3.5) 

where now stands for (0) for the sake of brevity. 

Let us evaluate {E(t)) by first integrating over while B^ is kept fixed. The first 
term in (3.5) is independent of and the second term is quadratic, leading to 



{E{t)) B = -f Tr[5„ B^f + ^Tr e2*A (A^ - d,d,) ^ (3.6) 



where only the leading order propagator is taken into account from the action. The integral 
and trace in the second term is given by, 

3(Ar2_i)^e-^'"'s*/^' = (3.7) 



using equation (B.2) from the appendix and where the ratio c = y/8t/L was introduced. 
The factor 3 comes from the trace over the Euclidean indices and the factor N'^ — 1 comes 
from the gauge trace. 

Let us now integrate over using the effective action (2.19). One needs to keep the 
tree level part only, all further terms are higher order in g^i. The second term in (3.6) is 
independent of B^ while for the first term we need the matrix integral 

12 ' 



/ dB^TrlB^, B,Y exp (^Tr[S,„ B,\ 



Mtvtr R 12^ 



(3.8) 



/dBexp(^^TY[S^,S,] 

Even though the integral is quartic it can easily be done with the result A^^ — 1 essentially 
determined by the dimensionality of the integral. Combining (3.6), (3.7) and (3.8) we 
obtain. 
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Figure 1. Finite volume correction factor S(c). 

which is the advertised final result (1.2). The finite volume correction term 6{c) is plotted 
on figure 1 as a function of the ratio c. As can be seen the correction never reaches 10% 
for < c < 1/2. 



4 Running coupling 

The result (3.9) can be used to define a non-perturbative running coupling scheme in which 
the running scale is = 1/L. As one changes the scale one keeps c fixed. Then the scheme 
is defined by the couphng constant 



3{N^ - + 6{c)) 



(4.1) 



where now the expectation value on the right hand side is understood non-perturbatively. 
The results from the preceding sections ensure that the above defined coupling for small L 
will run according to the universal 1-loop /3-function. Different choices for c correspond to 
different schemes. 

A note is in order about the 2-loop /3-function. As is well known both the 1 and 2 
loop coefficients are universal under a scheme change of the type g = g{l + 0{g^)) where 
the expansion on the right hand side only contains even powers of the coupling. However 
if one allows scheme changes of the type g = g{l + 0{g)) where the expansion contains 
both even and odd powers then only the 1-loop coefficient remains scheme independent. 
Our scheme is related to the MS scheme by such an expansion since it is easy to see that 
both even and odd powers of the coupling will appear as subleading terms to the leading 
result (1.2) but fractional powers will not. Our scheme is nevertheless well-defined and has 
for instance the property that if a theory has an infrared fixed point in one scheme it will 
have a fixed point in our scheme as well. 
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Figure 2. Discrete /^-function of SU{3) gauge tlieory coupled to Nf — 4 flavors of massless 
fundamental fermions for a scale change of s = 3/2. The results at 3 lattice spacings are shown 
together with the continuum 1 and 2-loop result from (5.3) for comparison. 

In order for the system to be controlled by a single scale L the bare fermion mass was 
set to zero in the preceding sections. The spectrum of the Dirac operator nevertheless has 
a gap ~ 1/L due to the non-trivial boundary conditions for the fermions. 

5 Numerical results 

We have tested the new running coupling scheme in SU (3) gauge theory coupled to Nf = 4 
massless fundamental fermions. The Schrodinger functional analysis of the same model can 
be found in [19, 20]. The fermion action was the 4-step stout improved [21] staggered action 
with smearing parameter g = 0.12. Since the number of flavors is a multiple of four no 
rooting was necessary. For the gauge sector tree level improved Symanzik action [22, 23] 
was used. The hybrid Monte Carlo algorithm [24] was used together with multiple time 
scales [25] and Omelyan integrator [26]. 

The observable E{t) and the flow itself can be discretized in a number of ways. Both 
the discretization in [3] and also the tree level improved Symanzik discretization of [7] 
was measured. We have found that the latter displays better scaling as expected hence 
in the following only the results from the Symanzik discretization will be presented. The 
bare quark mass was set to zero and anti-periodic boundary conditions were used for the 
fermions in all four directions. As mentioned in the previous section this leads to a gap 
~ 1/L in the spectrum of the Dirac operator. The gauge field was periodic in all directions. 
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2.58(2) 


1.77(1) 


1.352(8) 



Table 1. Measured renormalized couplings gl{L) from (4.1) at c = 0.3 and given bare couplings /3 
and lattice volumes L/a. 



The choice of0<c<l/2is limited by the observations that a small c leads to 
large cut-off effects while large c leads to large statistical errors. We found that c = 0.3 
is a convenient choice and from here on will drop the index c or i? on the renormalized 
coupling g^. 

The discrete version of the /3-function, or step scaling function, was computed for a 
scale change of s = 3/2. Three lattice spacings are used corresponding to 12^ — t- 18^, 
164 ^ 24^ and 24"^ 36^. Then the discrete /3-function 



g^{sL)-g^{L) 

log(s2) 



(5.1) 



can be calculated as a function of g'^{L). Holding L fixed in physical units the continuum 
limit corresponds to L/a ^ oo. 

The numerical results can be compared with the perturbative /3-function for small 
renormalized couplings. The 2- loop /3-function is given by 



L- 



:dg' 
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9 



+ b2 
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h 25 



154 
~3~ 



dL2 n67r2 ^(167r2)2 ' 
The discrete /3-function up to 2 loops for a finite scale change s is then 



g\sL)-g^{L) 

log(s2) 



^^16^ + 



9\L) 
(167r2)2 



(5.2) 



(5.3) 



which will be used for comparison although the zero mode of our finite volume scheme will 
introduce modifications which have not yet been calculated. 

The measured results for the renormalized coupling at each bare coupling and lattice 
volume are tabulated in table 1. At the volumes 12"^, 16^^, 18^, 24^ and 36^ the number of 
equilibrium trajectories were 10000, 10000, 10000, 8000 and 4000, respectively and every 
10*^ configuration was used for measurements. Auto correlation times were also measured 
and are around 10 — 30, 10 — 40, 10 — 70, 30 — 100, 30— 100 for the five volumes, respectively. 
The lower auto correlation times in the indicated intervals correspond to larger /3 and the 
higher ones to smaller /3. 

The discrete /3-function obtained from the data is shown on figure 2. The continuum 
extrapolation can be performed in (at least) two different ways. In the first method a 
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Figure 3. Paranictrization of the curves g'^{f3) at fixed lattice volumes using the expression (5.4). 
Red: 12*, green: 16^, dark blue: 18*, magenta: 24*, light blue: 36". 

cubic spline interpolation is done at fixed L/a — )• sL/a for {g^{sL) — {L)) / \og{s^) as a 
function of g'^{L). Then the resulting three curves together with their errors are used for 
the continuum limit at each fixed g'^{L). The continuum extrapolation is linear in jlT' 
since both the action and the observable only contain O(a^) corrections. This latter step 
is repeated for each value of g^{L). 

In the second method, similarly to [19], the dependence of g'^{/3) on f3 at fixed L/a is 
parametrized by the expression 



and the coefficients Cm are fixed by fitting to the measured values. The x^/dof values from 
the fits for the five volumes are 1.59, 0.39, 0.45, 1.11 and 0.08, respectively from 12"^ to 36''. 
The fitted curves together with the data are shown on figure 3. Since the parametrization 
is linear in the coefficients Cm the error on the fitted curve can be computed in a straightfor- 
ward manner. Then g'^{L) together with the discrete /3-function {g'^{sL) — g^ (L)) / log{s'^) 
and its error can be obtained for any /3 for all three lattice spacings corresponding to 
12^ — )• 18^, 16^ — )• 24^^ and 24^ — )• 36^. From here the procedure is identical to the pre- 
vious method; at fixed g'^{L) the three discrete /3-function values are extrapolated to the 
continuum assuming 0{a^ /Lp') corrections. 

The continuum extrapolation is shown on figure 4 for both methods and for four 
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Figure 4. Continuum extrapolations of tlie discrete /?- function for four selected g^{L) values 1.4, 
2.2, 3.0 and 3.8. Both methods are shown together with the /dof values of the fits. 

representative values of g'^{L), 1.4, 2.2, 3.0 and 3.8 together with the /dof values of the 
fits. The continuum results agree nicely between the two methods. 

It is reassuring to note that the continuum extrapolations from the two methods yield 
continuum results that agree with each other within error showing the robustness of the 
procedures. Also the continuum result is quite insensitive to the order of the polynomial 
used in (5.4) or other details of the fitting procedures. 

The final continuum extrapolated result agrees approximately with the 2-loop pertur- 
bative expression (5.3) as shown on figure 5 (only the final result from the first method 
is shown, but the second one gives a result which agrees with it within errors in the en- 
tire g^{L) range). As noted in section 4 our scheme is related to the MS scheme via 
gl = ^1^2(1 + fli(c)ff^jg + . . .) where ai(c) is non-zero leaving only the first /3-function co- 
efficient scheme independent. It can be shown from the measured gradient flow at c = 0.2 
that the discrete /3-function in figure 5 is not sensitive to the volume beyond the leading 
(5(c) correction factor. This explains the approximate agreement with the 2-loop universal 
/3-function keeping contributions from ai(c) undetectable within errors. 



- 12 - 




Figure 5. Discrete /3-function of SU{3) gauge tlieory coupled to Nj ~ 4 flavors of massless 
fundamental ferniions for a scale change of s = 3/2. The continuum extrapolated result from 
method 1 (see text for details) is shown together with the 1 and 2-loop results from (5.3) for 
comparison. 

6 Conclusion and outlook 

The Yang-Mills gradient flow - or as implemented on the lattice, the Wilson or Symanzik 
flow - is a promising tool for lattice gauge theory. In order to use it for a running coupling 
scheme where the running scale is the size of the finite 4-dimensional box one needs to 
compute the flow in perturbation theory at finite volume. In this work the necessary 
calculations were done in dimensional regularization and the MS scheme. The fact that 
the volume is finite necessitates the separation of the gauge Fourier modes into zero and 
non-zero modes. The non-zero modes can be treated in 1-loop perturbation theory while 
the non-trivially interacting zero modes need to be treated exactly. The result at leading 
order contains both algebraically and exponentially suppressed finite volume correction 
terms relative to the infinite volume result. 

The flow was then used to define a new scheme for the renormalized running coupling, 
agreeing with all other schemes at small couplings, as it should. The new non-perturbatively 
well-defined scheme is actually a one parameter family of couplings, all behaving universally 
for small values. The free parameter can be optimized for various targets such as small 
cut-off effects and/or small statistical errors. Further advantages of the scheme is the 
fact that the necessary observable can be evaluated in a Monte Carlo simulation at small 
computational cost relative to the HMC algorithm itself because only gluonic observables 
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are needed. In addition no extrapolation is needed for large Euclidean times. 

The flow in infinite volume, as noted in [3], can also be used for a running coupling def- 
inition. In this setup the coupling runs with the scale /j, = Lattice implementation 
of this running over many orders of magnitudes requires additional control to keep finite 
volume effects small. The scheme presented in this work circumvents this problem because 
the running scale is the volume itself, similarly to the Schrodinger functional method [9]. 

It would be very interesting to calculate further terms in the expansion (1.2) as well 
as the leading cut-off effects to it. Also, possible tunneling events at small renormalized 
coupling need to be investigated in the future but in the numerical simulations so far we 
have not observed any for the measured observable. 
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A Conventions 

The gauge field is taken to be anti-hermitian and the covariant derivative and field strength 
tensor are given by 

Ff,u = d^A^ - d^Af, + [A^, A^] (A.l) 

The SU{N) trace is then negative definite. The notation D^{B) = + B^^^ is also used 
and if it acts on gauge fields it is understood to be in the adjoint representation. When 
acting on fermions it acts in the representation R. 

The Dirac operator ij) = 'jfj.D^ is also implicitly acting in representation R. For the 
commutators of gamma matrices a^j/ = ^{'j/^'jiy — Jul/i) is used. 

B Evaluation of momentum sums 

The Jacobi elliptic function is used in the text, 

oo 
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which has the property 

1} (e"'^/*) = Vi^ (e-^*) . (B.2) 
Let us introduce h{t) = i?(exp(— ivr)). Then we have, using property (B.2), 



dt + (^h^^t) - 1)) > (B.3) 



tt' ( 2d 



r(s) V2s(2s - d) 

which leads to the finite value 5 = —5.545177 for s = 1 and d = 4. For the other two 
constants we similarly have, 



„4 „s/-Z rco 



and using again property (B.2) we obtain for s = 4 and d = 4 — 2e, 

2 2 

a = ^ + 1.622255 + 0(e), /3 = ^ - 0.063112 + 0(e) . (B.5) 
8e 246: 

Clearly, a — 3/3 is finite. The last remaining constant 7 is not independent from the rest, 
we have 7 = da + d{d — l)/3, leading to 

7 = — - 2.492991 + 0(e) . (B.6) 
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